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In a leading order pQCD model, we have studied the heavy flavor production in p+p collisions at 
RHIC and LHC energy. Leading order pQCD models require a K-factor. At RHIC energy, ^/s—200 
GeV, we fix K such that the model reproduces the integrated charm yield, dN cc /dy, estimated by 
the STAR and the PHENIX collaboration in p+p collisions. The model then explains the STAR 
data on the transverse momentum distribution of open charm mesons (D°) and decay electrons in 
p+p collisions. The p+p predictions, scaled by the number of binary collisions, also explain the 
electron spectra in STAR p+d collisions and PHENIX Au+Au collisions in different centrality bins. 

iy-} ■ Assuming that at LHC energy K-factor is of the order of unity, we have used the model to predict 

the transverse momentum distribution of D and B mesons and also of electrons from semileptonic 

f^) ■ decay of D — > e and B — » e, in p+p collisions at LHC energy, y^i=14 TeV. 
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I. INTRODUCTION 

Recently STAR collaboration has measured the transverse momentum distribution of open charm mesons, D°(D°) 
from direct reconstruction of D°(D t) ) — > K^tt in d+Au collisions and indirect electron-position measurements via 
charm semileptonic decays in p+p and d+Au collisions at -\/s=200 GeV [1]. Electron spectra show approximate binary 
collision scaling between p+p and d+Au collisions. The total charm cross section per nucleon-nucleon interaction for 
d+Au collisions at y / s=200 GeV is a cc = 1.3+0.2+0.4 mb, considerably larger than the standard PYTHIA prediction. 
Semileptonic decay spectra of electrons has also been measured by the PHENIX collaboration [2] in p+p and in 
Au+Au collisions, at ^—200 GeV. In Au+Au collisions, PHENIX measured charm electrons in different centrality 
bins. PHENIX data also indicate that the binary collision scaling holds for charm production in Au+Au collisions. 
However, PHENIX measured lesser number of charm. Their estimate, a cc = 622 ± 57 ± 160 ^b, is approximately 
one half of the STAR estimate. PHENIX Au+Au charm results are in contrast with the strong violation of binary 
collision scaling observed in non-charmed (light) hadrons. All the four RHIC experiments, STAR, PHENIX, PHOBOS 
and BRAHMS reported large suppression of high pt (light) hadrons in Au+Au collisions [3-6]. High pt suppression 
of non-charmed (Light) hadrons are generally explained in terms of partonic energy loss in a dense medium [7,8]. 
Indeed, they are the most important inputs for the claim that very high density medium is created at RHIC Au+Au 
collisions. However the charm results are not contradictory with the high px suppression observed in non-charmed 
hadrons. Energy loss do not reduce the number of partons, rather change the px distribution. Thus if in Au+Au 
collisions, heavy quarks suffered energy loss, momentum distribution of charmed electrons would not have scaled with 
charmed electron distribution in pp collisions. The PHENIX Au+Au data on charm electrons then confirm the believe 
that heavy quarks suffer no or little energy loss in the medium. 

Theoretical understanding of heavy flavor production is of considerable interest. Perturbative QCD is better appli- 
cable for heavy flavor production. Their large mass provide a natural scale for perturbative expansion. Experimental 
data on heavy flavor production can test the perturbative models and provide important inputs like various mixing 
angles. In heavy ion collisions at RHIC, heavy flavor production can give important information about the initial 
condition of the medium produced. Due to their large mass, heavy flavors are produced in initial hard collisions 
and in a short time scale and are ideal probe for the initial condition of the medium produced. Passage of a heavy 
flavor through a deconfmed medium can alter the momentum distribution of heavy flavored hadrons. Semileptonic 
decay of a heavy flavored hadron D(B) — > Xev, add to the charmonium background, an important signal of the 
confinement-deconfinement phase transition. Open charm yield is also important for understanding charmonium pro- 
duction. Quite early, Matsui and Satz [9] suggested that charmoniums will be suppressed in a deconfmed medium. 
However, recently there have been suggestions that rather than suppression, at RHIC, charmonium yield will be en- 
hanced due to recombination effect [10]. Open charm yield can help to decipher the issue of suppression/enhancement 
of charmoniums. 

Heavy flavor production in pQCD models has been studied earlier also [11-13]. Recently Cacciari et al [13] made 
an up-to-date prediction for open charm and bottom production at RHIC p+p collisions, using the First Order plus 
Next-to-leading-Log (FONLL) level. In the present paper, we have studied heavy flavor production in a leading order 
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pQCD model. Leading order pQCD models are simpler but need a if-factor. K factor accounts for the neglect of all 
the higher order terms in the perturbative expansion. K-factor depends on energy. Energy dependence of if-factor 
has been studied in [14]. i-C-factor decreases with energy and at RHIC energy, K ~ 3.4. At still higher energy K 
factor approaches unity. In the present paper, at RHIC energy K is obtained by fitting the integrated charm yield 
dN c£ /dy in STAR and PHENIX p+p collisions. With the K-factor so fixed, the model is used to explain the STAR and 
PHENIX data on the transverse momentum distribution of D-mesons and decay electrons. As it will be shown later, 
the model, with heavy quark fragmentation function, parameterised by Braaten et al [15], give excellent description 
to STAR and PHENIX data. The more well known Peterson's fragmentation function [16] undcrpredict the data by 
a factor of 3 or so. The model can serve as a baseline for comparing heavy flavor data in RHIC Au+Au collisions and 
in future LHC energy collisions. 

The paper is organised as follows: in section 2, we briefly describe the pQCD model. Results for heavy flavor 
production and comparison with RHIC experiments (STAR and PHENIX) are done in section 3. Model predictions 
for heavy meson production and decay electrons at LHC energy, v / s=14 TeV p+p collisions are also given in sec. 3. 
Summary and conclusions are given in section 4. 



II. PCQD MODEL OF CHARM QUARK AND CHARMED MESON PRODUCTION 

In pQCD models, semileptonic decay of a heavy meson, occur in three step: (i) production of a heavy quark in 
parton level collisions, (ii) fragmentation of the heavy quark in to a heavy mesons and (iii) semileptonic decay of 
heavy meson. Schematically transverse momentum spectra of decay electrons in p+p collisions can be written as 
(schematically), 

where the symbol (g) denote generic convolution and D(Q — > Hq) and J(Hq — > e) represent the fragmentation of 
heavy quark Q into a heavy flavored hadron Hq, and semileptonic decay of Hq into electron respectively. 

Detail of the evaluation of Eq.l is given in the appendix. For evaluation of Eq.l, parton distribution functions and 
fragmentation function of a heavy quark to fragment into a heavy hadron, are required. For the parton distributions, we 
have used CTEQ5L parameterisation. Several fragmentation function have been used in literature [15-18]. Most well 
known is the parameterisation by Peterson et al [16]. Peterson's parameterisation does not contain spin information 
and is parameterised as, 

N 

DC/C[Z) = z[l-l/z- eQ /(l-*)]' (2) 

here N is the normalisation and eq is a parameter. The parameter sq is approximately the ratio of square of 
constituent light quark (q) and heavy quark (Q) masses, sq = m 2 /mQ. eq — m^/mg depend on the heavy quark 
mass. Most commonly used values are, eq = 0.5 for fragmentation of a charm quark into a D-meson and eq = 0.006 
for fragmentation of a bottom quark into a B-meson. More recently Braaten et al [15] parameterized the heavy 
quark fragmentation function. It contains the spin information. The pseudo-scalar and vector meson fragmentation 
functions arc written as, 



Dq^p(z) = N- y '—^[6 - 18(1 - 2r)z + (21 - 74r + 68r 2 )z 2 - 2(1 - r)(6 - 19r + 18r 2 )^ + 

[1 — (1 — r)z\ 

3(l-r) 2 (l-2r + 2r 2 )z 4 ] (3) 
rzi\ - z) 2 

Dq^ v (z) = N- ( - '-—[2 - 2(3 - 2r)z + 3(3 - 2r + Ar 2 )z 2 - 2(1 - r)(4 - r + 2r 2 )z 3 + 

[1 — (1 — r)z\° 

(l-r) 2 (3-2r + 2r 2 )z 4 ] (4) 

Parameter r is approximately the ratio of constituent quarks masses, r = m q /m,Q. r = 0.1 for c-quark fragmentation 
to D-mesons, and r = 0.03 for fragmentation of 6-quarks into a B-meson. In the present work we have used the pQCD 
pseudo-scalar fragmentation functions parameterised by Braaten et al [15], henceforth call FF-I. For comparison, we 
will present results obtained with the fragmentation function parameterised by Peterson et al [16] (henceforth call 
FF-II). As will be shown later, fragmentation function FF-II yields are nearly a factor of two less than the yield 
obtained from FF-I. Peterson's fragmentation function is comparatively hard. 
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Fragmentation functions, Eq.2 and 3 are used as the initial distribution at momentum scale, q 2 = Tn 2 , b y In 
the present work, we have used m c =1.5 GeV and mf,=4.4 GeV as charm and bottom quark mass respectively. 
Fragmentation function at other momentum scales are obtained from numerically integrating the Alterelli- Parisi 
evolution equation, 

d f du z 

V 2 tt^D Qh (z,v 2 )= -^P Q ^ Q (-,fi)D Q ^ H (y,^) (5) 

<V Jz y v 

where Pq^q is the splitting function: 

P Q ^ Q {z,f) = ^±(——A (6) 
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FIG. 1. (a) Solid line is the pQCD prediction for 
the pt distribution of c-quarks at y / s=200 GeV with 
K=l. The dashed and short-dashed lines are the 
model predictions for _D°-mesons and decay electrons 
with FF-I. Red dashed line and short dashed lines are 
the model predictions for D-meson and decay elec- 
trons with FF-II. (b) same for b-quarks, B-mesons and 
B — > e decay. 
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D->e (FF-I) 

B->e(FF-l) 

total(FF-l) 

D->e(FF-ll) 

- B->e(FF^II) 

- total(FF-M) 




FIG. 2. PHENIX measurements of pr distribution 
of charm electrons in p+p collisions at RHIC energy, 
^=200 GeV. The pQCD predictions with FF-I , are 
shown as dash-dot and dash-dot-dot lines, for D and 
B meson decay respectively. The solid line is the sum 
of the two contributions. Results obtained with FF-II 
are shown as red lines. The K-factor is 1.8. 



III. RESULTS 



A. Transverse momentum distribution of D ,B and semileptonic decay electrons at RHIC 



The pQCD model predictions for the transverse momentum distribution of heavy mesons and decay electrons in 
y / s=200 GeV p+p collisions are shown in Fig. la and lb. Here we have taken K=l. In Fig. la, solid black line is 
the pt spectrum for charm quarks. Charm quark fragmentation to D mesons and semileptonic decay of D-mesons, 
obtained with the fragmentation function FF-I is shown as black dashed and short-dashed lines respectively. The same 
obtained with the fragmentation function FF-II are shown as red dashed and red short-dashed lines, px distribution 
of D-mesons closely scale with c-quarks but pt distribution of electrons falls much faster than that of c-quarks 
or of D-mesons. We also note that yield of D mesons and also of decay electrons, obtained with the fragmentation 
function FF-I are approximately a factor of 2 larger than the yield obtained with the parameterisation FF-II. Peterson 
fragmentation function (FF-II) is comparatively hard. Same results, for bottom production are shown in Fig. lb. At 
RHIC energy, low pt heavy quarks are dominantly charm. Heavy mass of bottom quarks inhibit their production 
at low pt- However, at large pt, bottom production evenly competes with charm production. Consequently, at 
large pt semileptonic decay electrons will have contributions from D as well as from B mesons. For B-mesons also 
fragmentation function FF-II produces factor of two less number of B-mesons than the fragmentation function FF-I. 
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As mentioned earlier, leading order pQCD models require a k-factor. STAR and PHENIX collaboration had 
estimated the total charm yield in p+p collisions at \/s=200 GeV. While the PHENIX collaboration [2] obtained 
da c5 /dy = 143 ± 13 ± 36 fib in the p T range, OAGeV < p T < 4GeV, the STAR collaboration [1] obtained larger 
value, da c5 /dy = 0.29 ± 0.04 ± 0.08 mb in p T range, lGeV < p T < 4GeV. We have fixed the K factor such that the 
integrated charm yield in the model reproduces the central value of the PHENIX and STAR estimate. We obtain, 
K=1.8 and K=7.0 for PHENIX and STAR experiment respectively. We could have fixed K to reproduce the D meson 
or the electron yield. But fixing the K from charm yield has the advantage that the model can now test different 
fragmentation functions by comparing the model predictions with experiment. Before we continue we note that the 
large uncertainty in K (K=1.8 for PHENIX experimental data and 7.0 for STAR experimental data)is a direct result 
of large difference in integrated charm cross sections in two measurements. We hope that in future the large difference 
in charm cross sections, in two experiments, will be removed, and K can be fixed more accurately. 

With the only parameter of the model fixed, we can compare model predictions with other experimental observable, 
namely, pt distribution of D° and decay electrons in STAR and PHENIX experiments. In Fig. 2, pt distribution of 
charm electrons in pp collisions, measured by the PHENIX collaboration [4] at -^=200 GeV is shown. Present pQCD 
model predictions obtained with the fragmentation functions FF-I (black lines) and FF-II (red lines) are also shown 
in Fig. 2. For both the fragmentation functions, the K-factor is fixed, K=1.8. We have shown the contribution of both 
D (the dash-dot line) and B mesons (dash-dot-dot line) decay. Solid lines are the sum of the two contributions. At 
low pt, decay of B-mesons contribute negligibly to the total electron spectrum. For pt <2 GeV, bottom contribute 
less than 10% to the total yield. In the pr range 2-3 GeV bottom contribute 10-20% and in px range 3-4 GeV bottom 
contribution rises to 20-30%. This numbers compare favorably with PHYTHIA estimate [19]. Present pQCD model 
with FF-I well reproduces the PHENIX data. On the other hand fragmentation function FF-II, though reproduce the 
shape of the spectra, underpredict the PHENIX data by a factor of ~2. This is understood. As mentioned earlier, 
Peterson's fragmentation functions are hard and it produces less number of heavy mesons. 
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FIG. 3. Transverse momentum distribution of 
charmed electrons in Au+Au collisions at RHIC en- 
ergy, ^5=200 GeV. Solid lines are pQCD predictions 
(K=1.8) for p+p collisions, scaled by the nuclear over- 
lap function. 
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FIG. 4. STAR measurements of pr distribution 
D-meson (solid circles) and charmed electrons (solid 
squares) in d+Au collisions at RHIC energy y/s—200 
GeV. STAR pp results for decay electrons are shown 
as solid triangle. The solid and dashed lines are the 
pQCD predictions with FF-I and FF-II respectively 
with K=7.0 . 



In Fig. 3, PHENIX measurements [2] in Au+Au collisions, in different ccntrality bins are shown. As mentioned 
earlier, charm electrons scale with binary collision numbers. In Fig. 3, present pQCD model results, with K=1.8 are 
shown. The solid and dashed lines are obtained with fragmentation function FF-I and FF-II respectively. We have 
included the B meson decay contributions. For comparing with PHENIX data the pp-predictions are multiplied by the 
nuclear overlap function. As expected, model predictions with Braaten parameterisation of fragmentation function 
agree well with the experiment. Here again, Peterson's parameterisation underpredict the data approximately by 
a factor of 2. Analysis of PHENIX data confirms that heavy quarks suffer little or no energy loss in the medium 
produced in Au+Au collisions at RHIC energy. 
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We now compare the present model predictions with the STAR experiment. STAR collaboration [1] have measured 
PT-distribution of D-mesons in d+Au collisions at GeV. They have also measured the semileptonic decay 

electrons in d+Au and p+p collisions. In Fig. 4, STAR data are shown. As mentioned earlier, for STAR experiment, 
we use K—7.0 . In Fig. 4, solid and dashed lines present model predictions with fragmentation functions FF-I and 
FF-II respectively. For d+Au collisions, we have multiplied the results for pp predictions by N^n = 7.5. We note that 
the fragmentation function FF-I, with K—7 can describe the STAR data on charm production in p+p, p+d collisions. 
FF-II, as before leads to poorer description of the data. 



B. Transverse momentum distribution of D ,B and semileptonic decay electrons at LHC 



We have shown that the present pQCD model with the heavy quark fragmentation function, parameterised by 
Braaten et al [15], reproduces the STAR and PHENIX data on the transverse momentum distribution of D meson 
and semileptonic decay electrons in y / s=200 GeV p+p, d+Au and Au+Au collisions. With a reasonable estimate of 
K, the model can be used to predict heavy flavor production in pp collisions at LHC energy, As mentioned earlier, K 
factor decreases with energy and at LHC energy, expected value is K ~ 1 — 1.5 [14]. In the present section, we give 
predictions for the transverse momentum distribution of D, B and semileptonic decay electrons with K=l. We only 
show the predictions obtained with the fragmentation function FF-I, parameterised by Braaten et al [15]. As shown 
earlier, RHIC data are better explained with this parameterisation. Before we proceed, a note of caution is in order. 
The model is tested against RHIC data, limited to pr < 4.5 GeV. Accuracy of the model beyond pr=4.5 GeV is not 
tested. In the following, we will show model prediction over a wide range of pr, but the predictions for pr > 5 GeV 
may not be accurate. 
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FIG. 5. Transverse momentum distribution of D 
mesons (solid line) and B mesons (dashed line) in pp 
collision at LHC energy, ^5=14 TeV. 
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FIG. 6. Transverse momentum distribution of elec- 
trons in semileptonic decay of D mesons (dash-dot 
line), B mesons (dash-dot-dot line) in pp collisions at 
LHC energy ^=14 TeV. The solid line is the sum of 
the two contributions. 



In Fig. 5, model predictions, for the transverse momentum distributions of D and B mesons in pp collisions at LHC 
energy y / s=14 TeV, are shown. As it was at RHIC energy, at LHC energy also, low px heavy mesons are dominated 
by the D-mesons. At low pr B-meson yield is significantly lower than the D meson yield. Only at large pr, B- meson 
contribute significantly to total heavy meson production. However, compared to RHIC energy (see Fig.l), at LHC 
energy both D and B meson yields are increased. High pt yields are significantly increased, by a factor, 1000 or more. 

In Fig. 6, pt distribution of electrons from semileptonic decay of D and B meson are shown. Low pt electron spectra 
is dominated by the electrons from decay of D-meson. Below pr=l GeV, B mesons contribute less than 10% to the 
total yield. B meson contribution increases rapidly with pr and in pr range 2-4 GeV, it rises to 30-50%. However, at 
very large pr > 6GeV B mesons contribution decreases again. 



5 



IV. SUMMARY AND CONCLUSIONS 



In a pQCD based model, we have studied the heavy flavor production at RHIC energy pp collisions. At RHIC 
energy, heavy flavor production is dominated by the charm quarks. Below pt~2 GeV, more than 90% heavy quarks 
are charm. Production of bottom quark picks up only at large pr and at pt—5 GeV, it contribute to 20% of to 
total heavy quark yield. As a consequence of dominance of charm quark at low px, at RHIC low px heavy mesons 
are dominantly D-mesons. Similarly decay electrons are dominantly from decay of D meson. We have compared the 
pQCD model predictions for D and decay electrons with STAR and PHENIX data. This models require a K-factor. At 
RHIC energy, we fix the K-factor to reproduce the integrated charm yield d<r cc /dy, estimated in STAR and PHENIX 
experiments. With K-factor fixed from total charm yield, the model, with fragmentation function parameterised 
by Braaten et al [15], reproduces the transverse momentum spectra of semileptonic decay electrons measured by 
the PHENIX and STAR collaboration in ^=200 GeV p+p collisions, pp predictions scaled by the nuclear overlap 
function also explain the PHENIX Au+Au charm data in different centrality bins, pt distribution of D° in p+d 
collisions, measured by the STAR collaboration is also explained in the model. The results confirms that charm 
production in d+Au and Au+Au collisions obey binary collision scaling law. The other, much used fragmentation 
function parameterised by Peterson et al [16], though reproduces the shape of the distributions, undcrprcdict the 
yields by a factor of ~ 2-3. Assuming that at LHC energy, y / s=14 TeV, K-factor approached unity, we have used the 
model to predict heavy flavor production in p+p collisions at LHC energy. Compare to RHIC energy, at LHC energy 
collisions, heavy flavor production is increased. At large px increase is more than a factor 1000. 



APPENDIX A: PQCD MODEL FOR P T DISTRIBUTION OF SEMILEPTONIC DECAY OF D — > XEv 

Production cross section of a heavy quark c (charm) in a pp collision can be written as [20] 

partons „ . 

da v~~* / , , i. / ^7, „ , s da 



dyd?pj 



K dx a dx b f i (x a ,Q 2 )f j (x b ,Q 2 )-^(ab^cd)5(s + i+u-ml-ml-ml-mf) (Al) 

Z — ' J 7T dt 



where x a and Xb are the fractional momenta of the colliding partons. (s,t, u) are the Mandelstam variables for the 

da 

di 



subprocess ab — > cd. % is the subprocess cross-section. The factor K takes into account the neglect of higher order 



terms. 

For a open charm production subprocesses ab — > cd of interest are: 

99 -> CC (A2) 

qq -» CC (A3) 

gQ -> gC (A4) 

qQ -» qC (A5) 

where C denote a charm quark, g and q denote a gluon and a light quark respectively. 

Matrix elements for these processes are calculated in [21]. For charm quark mass M subprocess cross sections are, 



, ,, n2u 12,„, - w ,, 2 „ x 8(M 2 -t)(M 2 -u)-2M 2 (M 2 +t) , . , 

W\ 2 gg ^= ^ 2 (Q 2 )[^(M 2 - t){M 2 -u) + - K - (M 2_£)2 " ( A6b ) 

8{M 2 -t){M 2 -u)-2M 2 (M 2 +u) 2 M 2 {s-AM 2 ) 



+ 



3 (M 2 - u) 2 3 (M 2 - t) (M 2 - u) 

(M 2 - t)(M 2 -u) + M 2 {u - t) _ (M 2 - t) (M 2 -u) + M 2 {t-u) 
s(M 2 - 1) s(M 2 - it) J 

lA.fl 2 6 4 2 n ,2 rn2 ^ (M 2 u) 2 + (j - M 2 ) 2 + 2M 2 t 

\ M \qc^qc= Y n a W ' {2 ( A6C ' ) 

lA^I 2 A 2 (n 2, M.s-M 2 ){M 2 -u) U (s-M 2 ){M 2 -u) + 2M 2 {s + M 2 ) 

\M\ gq ^ gc = 7T a (Q )[ p + — (s-M 2 Y ( ' 



G 



+ 64(s- M 2 )(M 2 - u) + 2M 2 (M 2 +u) | 16 M 2 {AM 2 -i) 



{M 2 -u) 2 9 {s-M 2 ){M 2 -u) 

(a - M 2 )(M 2 -u) + M 2 (s - fi) (g - M 2 )(M 2 - u) + M 2 (s - fl) 
+ t(s - M 2 ) 1 £(M 2 - u) J 

For the subprocesses gg — > CC and — > C(7, we assume the following forms for the parton momenta a,b and c. 

Pa = (—^-,0,0, ) (A7a) 

Pb = (^— ,0,0, ) (A7b) 



Vc = {M T chy,p T , 0, M T shy), M T = ^M 2 +p 2 T (A7c) 

Mandelstam variables s = (p a +Pfc) 2 , t = (p ~ Pc) 2 and u = (p& — p c ) 2 , for these subprocesses can be calculated as, 

s = x a XbS (A8a) 

f = M 2 - x a y/sM T exp(-y) (A8b) 

tt = M 2 - x b VsM T exp(y) (A8c) 

For the subprocesses involving heavy quark excitations, we assume the following forms for p a ,Pb and p c 

Pa = (^,0,0,^) (A9a) 
z s 



P6 = (^ + ^.°.°.-^) ( A9b ) 

p c = (M T chy, p T , 0, M T s%) (A9c) 



Corresponding Mandelstam variables are 

s = Af 2 + a; a x b s (A10a) 
i = M 2 - x a sfsM T exp{-y) (AlOb) 
it = 2M 2 - x b \fsM T exp{y) (AlOc) 

One of the integration in Eq.Al can be eliminated with the condition s + i + u~ 2M 2 , and charm production cross 
section can be obtained as, 

da c Pa ^ ns f 1 , tl „ 2w , „ 2 2 x a x b + M 2 /s da 



dyd 2 p 7 

where 



V- / 7 , / ^2^,/ ^2n 2 x a x b + Mj/8 da., , A11 , 



_ - AM* /a 

Xb - 2x a -x T e-v (M2) 
x T e v - AM 2 Is 

< m = — =F- (A13) 

2 — ire y 

where M g = for the subprocesses, gg(qq) — ► CC and M g = M c for the subprocesses, gC — > gC and qC qC . 

If D D /c(z, ii 2 ) is the fragmentation function for quark C to fragment into £) meson h, transverse momentum 
distribution of D arc obtained as, 



D [dz„ , 2 . da c 

^ = y^/^'^^v (Au) 

irtons pi i fl o _i_ A/f2 / j 

E / o $^vc(z, m 2 ) jf ^ <**./<(»., Q 2 )fj(xb, Q 2 )~ Xa * x b a _ Zey S ^ ab ~> cd > ( A15 ) 



dy h d 2 q- 

partons „\ 
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For a given hadron momentum (qx) and rapidity (yh), corresponding quantities {pr,y) of the fragmenting parton 
are obtained from the relations, 



^M* + q 2 chy h = z^M 2 +p 2 T chy (Al6a) 
^M 2 + q 2 shy h = ^M 2 +p 2 T chy (A16b) 

Semileptonic decay of heavy mesons D — > Xev has been studied in detail [22,23]. We will not elaborate on the 
procedure. Transverse momentum distribution of electrons from the decay of the D- mesons can be calculated as [22], 

dN e f dN D 

-j— = / ^H(p T ,p' T )dp' T (A17) 
dp T J dp' T 

with 



H(pi 



> Pt) -J 2\p T \p T .p>J^ ] (M8) 



where / is the rest frame distribution of decay electrons. We have used the following rest frame spectrum for the 
decay D -> Xev [22], 

i c/r 

f{E e ) = -—=wg{E e ) (A19) 

where 

= E 2 (M 2 D -M x -2M D E e 
yv 1 M D -2E e y ' 

96 

W = (1 - 8m 2 + 8m 6 - to 8 - 24m 4 In m)M% 

m = M x /M D (A22) 
We assumed that for D-mesons decay, Mx = M^=0.497 GeV. 

Above formulation is given in terms of charm quarks and D mesons, They are also applicable for bottom quark and 
B-meson production with appropriate change of masses. For semileptonic decay of B-mesons, B — > Xev, Mx = Md- 

Decay electron spectra are appropriately weighted by the branching ratios. For D — > e branching ratio is 0.103 and 
that for decay of B — > e branching ratio is 0.109 [13]. We note that B — > D — > e channel also contribute to total 
electron spectrum, but we have ignored it, as second order contribution. 
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